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Abstract 

As neutron stars spin-down and contract, the deconfinement phase transi- 
tion can continue to occur, resulting in energy release (so-called deconfinement 
heating) in case of the first-order phase transition. The thermal evolution 
of neutron stars is investigated to combine phase transition and the related 
energy release self-consistently. We find that the appearance of deconfine- 
ment heating during spin-down result in not only the cooling delay but also 
the increase of surface temperature of stars. For stars characterized by in- 
termediate and weak magnetic field strength, a period of increasing surface 
temperature could exist. Especially, a sharp jump in surface temperature can 
be produced as soon as quark matter appears in the core of stars with a weak 
magnetic field. We think that this may serve as evidence for the existence of 
deconfinement quark matter. The results show that deconfinement heating 
facilitates the emergence of such characteristic signature during the thermal 
evolution process of neutron stars. 
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1. Introduction 

Fundamental properties of supranuclear matter in the cores of neutron 
stars, such as the chemical composition and the equation of state, are still 
poorly known. Simulations of the thermal evolution of neutron stars con- 
fronted with soft X-ray, extreme UV, and optical observations of thermal 
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photon flux emitted from their surface provide the most valuable informa- 
tion about the dense matter in the interior of these stars. 

As the interior density gradually increases of neutron stars, deconfine- 
ment phase transition continuously takes place inducing not only structural 
changes but also energy release in case of a first-order phase transition. The 
generation of energy increases the internal energy of the star which is called 
deconfinement heating (Haensel & Zdunik [H,Yu & Zheng (2|,Kang & Zheng 
(if). The thermal evolution of neutron stars is connected with their spin- 
down and the resulting changes in structure and chemical composition (from 
nucleon matter to deconfined quark matter) have been investigated in our 
work. We have investigated the thermal evolution of neutron stars with 
such a deconfinement phase transition(Kang & Zheng (if). The results show 
that deconfinement heating delays dramatically the cooling of neutron stars, 
which have a higher surface temperature compared with traditional cooling 
for the same age stars. 

Many efforts are devoted to explore the observational signal of a decon- 
finement phase transition which have been suggested in the form of char- 
acteristic changes of observables, such as the pulse timing(Glendenning et 
al. j3],Chubarian et al. ji[ ,Poghosyan et al.0]), brightness (Dar &DeRujula 
7|)and surface temperaturefSchaab et al.jsj], Blaschke et al.[9],Yuan & Zhang 
lOj] , Stejner et al. [ll|) of the pulsars during their evolution. The surface 
temperature changes when the quark matter appears in the cores of neutron 
stars. We explore the deconfinement signature by studying the changes of 
surface temperature in the thermal evolution process of neutron stars. 

In this paper, we reinvestigate the thermal evolution of neutron stars and 
look for a characteristic change in the surface temperature with deconfine- 
ment heating. The released energy can also be estimated as a function of 
the change rate of the deconfinement baryon number using the parameterized 
approach. The neutron stars containing quark matter are called hybrid stars. 
We take Glendenning's hybrid stars model (Glendenning (lij ) based on the 
perturbation theory developed by Hartlejl3| to study the rotational evolu- 
tion structure of stars. The Argonne V18 + 5v + UIX* model(APR)(Akmal 



et al. 14|) of hadronic matter and the MIT bag model of quark matter are 



used to construct the model of stars, but the medium effect of quark matter 
has been considered in quasi-particle descript ion (S chert le et al. 15| ) . 
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2. Deconfinement phase transition and neutron stars structure 



Early works on deconfinement phase transition, the Maxwell construc- 
tion, show a sharp transition taking place between the two charge-neutral 
hadron and the quark phases(Baym & Chin [la]). In the 1990s, Glenden- 
niiiK (Glendenning 0, H) 

pointed out that this assumption was too re- 



strictive. More generally, the transition can occur through the formation of 
a mixed phase of hadron matter and quark matter, with the total charge 
neutrality being achieved by a positively charged amount of hadronic matter 
and a negatively charged amount of quark matter. Following Glendenning's 
model, we use a standard two-phase description of the equation of state(EOS) 
through which the hadron and quark phases are modelled separately. The 
resulting EOS of the mixed phase is obtained by imposing Gibbs's conditions 
for phase equilibrium with the constraint that the baryon number and the 
electric charge of the system are conserved to the neutron star matter. 

The Gibbs condition for chemical and mechanical equilibrium at zero 
temperature between the two phases reads 

PHp{Pn,Pe) = PQp{Pn,Pe) (1) 

where pnp is the pressure of confined hadron phase and pqp is the pressure 
of deconfined quark phase. 

The conservation laws can be imposed by introducing the quark fraction 
X defined as x — Vq/V. Only two independent chemical potentials remain 
according to the corresponding two conserved charges of the /3-equilibrium 
system. The total baryon number density Pb is 

N B 

PB = ~^ = XPQP + (1 - X)PHP (2) 

the total electrical charge is 



= ^ = XQqp + (1 - x)Qhp (3) 
and the total energy density is 

e = y = X£qp + (1 - x)tHP (4) 

Using the Eqs.(l), (2), (3)and (4), we can obtain the EOS of mixed phase 
matter. In describing the hadronic part of the neutron star, we adopt the 
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APR model (Akmal et al. [14] ). For the EOS, it is based on the models 
for the nucleon interaction with the inclusion of a parameterized three-body 
force and relativistic boost corrections. We use the EOS of an effective mass 



bag-model for the quark matter part of the neutron star(Schertle et al. 15]). 

In Fig.l, we show the model EOS with deconfmement transition, which is 
the typical scheme of a first-order transition. The phase transition construc- 
tion in a two-component system leads to a continuous increase in energy 
per baryon in the mixed phase with increasing density. It is well known 
that hadron matter is the most stable phase at lower densities, and that 
quark matter is the most favorite phase at higher densities. Meanwhile, the 
mixed phase has the lowest energy at intermediate densities. We choose 
the parameters for quark matter EOS with s quark mass m s = 150MeV, 
coupling constant g = 3 and different bag constant B = 85MeVfm~ 3 , 
B = l08MeVfm~ 3 , B = 136MeV 'fm~ 3 respectively. 

With the EOS presented above, we are ready to study the structural 
evolution of the rotating neutron stars. In this paper, we apply Hartle's 
approach (Hartle [l3[) as in Kang & Zheng [3] to investigate the structure 
of the stars. By treating a rotating star as a perturbation on a non-rotating 
star and by expanding the metric of an axially symmetric rotating star in 
even powers of the angular velocity Q, we can obtain the structure of the 
rotating stars. 

The resulting gravitational masses of neutron stars for B = l08MeVfm~ 3 
are shown in Fig. 2, as a function of central baryon density for static stars as 
well as for stars rotating with the maximum rotation frequency v^. The solid 
almost horizontal connect configurations with the same total baryon number. 
In order to explore the increase in central density due to spin down, we cre- 
ated sequences of neutron star models. A Model in particular sequence has 
the same constant baryon number, increasing central density and decreasing 
angular velocity. 

3. Deconfmement heating 

There is deconfmement heating production due to spin-down in neutron 
stars due to the nuclear matter continuously converting into quark matter. 
The released energy had been estimated as a function the rate of change 
the deconfmement baryon number using the parameterized approach(Kang 
& Zheng [3|). Recently we studied the mechanism of energy release in de- 



tail (Kang et al. [18] ). Through studying a random process of infinitesimal 
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compression for the mixed phase region, we can calculate the energy release 
per baryon using the following formula 



Se - 5e = (-p- - -^—)5p B , (5) 

OpB 0p B 

where denotes the enthalpy change per baryon . 

The deconfmement heating is coupled with the rotation evolution of neu- 
tron stars. Combining the energy change with the evolutionary structure of 



neutron stars, we get the total heat luminosity (Kang et al. 18) 19) ) 



f de 

H dec = I -^v(t)p B dV (6) 

where v is the rotation frequency of the star. The spin-down of stars is due 
to the magnetic dipole radiation. The rotation frequency is given by 

v = -^L^Vsin 2 ^ (7) 
31c 6 

where / is the stellar moment of inertia, p = \BB? is the magnetic dipole 
moment, and 6 is the inclination angle between the magnetic and rotational 
axes. We now combine the energy release behavior with spin-down. Our 
recent work shows that exact calculations agree well with the earlier order 
of magnitude estimates and supports the previous parameterized approach. 
Because the change of rotating stars is sufficiently slow, the formula of heat 
luminosity can be replaced by the simple parameterized form(Kang & Zheng 

!), 

H dec = q^vit) (8) 
dv 

where q, about O.lMeV, is the mean value of energy release in the mixed 
phase and Nq represents the baryon number of quarks in the interior of the 
star. 



4. Signal of quark deconfinement and thermal evolution of neutron 
stars 

The cooling of neutron stars could take place via two channels - neutrino 
emission from the entire star and thermal emission of photons through the 
transport of heat from the internal layers to the surface. Neutrino emission 
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is generated in numerous reactions in the interior of neutron stars, e.g. as 
reviewed by Page et al.j2o[. For the calculation of cooling of the hadronic 
part of the neutron star, we use the main processes including the nucleon 
direct Urea (NDU) and the nucleon modified Urca(NMU) and the nucleon 
bremsstrahlung(NB). For the quark matter, we consider the main process: 
the quark direct Urea (QDU) processes on unpaired quarks, the quark modi- 
fied Urea (QMU) and the quark bremsstrahlung(QB). For pure neutron stars, 
the direct Urea reaction (the most efficient) is allowed only at very high den- 
sities because it is impossible to satisfy the conservation of momentum unless 
the proton fraction exceed the value where both the charge neutrality and 
the triangle inequality can be observed (Lattimer et al. [21] ) . However, for 
neutron stars that contains quark matter, this is not so because charge neu- 
trality does not have to be conserved locally for a mixed phase. Hence the 
NDU process is active in the mixed phase. 

The calculation of the evolution of the thermal energy of neutron stars 
(heating and cooling) is achieved by coupling their rotating structure and de- 
confinement phase transition. Fig. 3 displays the central density of different 
masses neutron stars for B = 108MeVfm~ 3 , as a function of their rotational 
frequency. In Fig.3, the dotted horizontal line indicates the deconfined quark 
matter produced and dashed horizontal line indicate the NDU processes trig- 
gered. In the interior of these stars, the emerging of deconfined quark matter 
accompanies the gradual energy release leads to rising of surface tempera- 
ture, and appearance of the NDU process can result in a rapid decrease of 
the temperature during the spin-down. For M = 1.5, 1.55, 1.6M Q neutron 
stars, the deconfinement phase transition occurs during the thermal evolu- 
tion process which may lead to appearance of a characteristic signature. We 
will now discuss in detail the quark deconfinement signal. 

We combine the equation of thermal balance with the rotating structure 
equations of the stars(Kang & Zheng j3|, Hartle [l3j]) and rewrite the energy 
equation in the approximation of an isothermal interior(Glen Sz Sutherland 

0) 

CviTuv)^ = -L?{Ti,v) - L~(T.,v) (9) 

C v (T u v) = c(r,T)(l - ^M)-V2 47rr 2 dr (1Q) 

Jo r 

L?(T h v)= j R{V \{r,T){l- 2 -^)-y^H^dr (11) 
Jo r 
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Where T s is the effective surface temperature, Tj(£) = T(r, t) is the redshifted 
internal temperature; T(r, t) is the local internal temperature of matter, 
and $(r) is the metric function(describing gravitational redshift)(Yakovlev & 
Haensel |23[). Furthermore, L£°(Tj, t>)and CV(Tj,t>) are the total redshifted 
neutrino luminosity and the total stellar heat capacity, respectively, which 
are functions of rotation frequency and temperature; c(r, T) is the heat ca- 
pacity per unit volume. = AnR 2 {v)aTf{l — R g /R) is the surface photon 
luminosity as detected by a distant observer(i? ff is the stellar gravitational 
radius). The effective surface temperature which is detected by a distant 
observer is T s °° = T s a/I — R g /R. T s is obtained from the internal temper- 



ature by assuming an envelope model(Gudmundsson et al. [24|, Potekhin 
et al |25|). Using to Eqs.(7)-(11), we can simulate the thermal evolution of 
neutron stars with deconfinement heating. 

In Fig.4 and Fig.5, we present thermal evolution behavior of a 1.6M Q 
neutron star for different magnetic fields (10 9 — 10 12 G). Due to the coupling 
of thermal evolution and spin-down, all curves(with deconfinement heating 
and without deconfinement heating) show clear magnetic field dependence. 
During the star's spin down, deconfined quark matter appears in the core of 
the star at a spin frequency of v — 1123Hz, the surface temperature drops 
rapidly when the NDU process(enhanced cooling) occurs at a spin frequency 
of v = 492Hz. It is evident that the temperature of the curves with de- 
confinement heating(solid curves) are higher than for the standard cooling 
scenario without deconfinement (dotted curves). We can observe a competi- 
tion between cooling and heating processes from the heating curves, where 
deconfinement heating can produce a characteristic rise of surface tempera- 
ture and even dominate the history of thermal evolution. Eventually, they 
reach a thermal equilibrium, where the heat generated is radiated away at 
the same rate from the star surface. We find the weaker magnetic fields have 
the larger change of temperature. The low magnetic field (10 9 G) produces a 
sharp jump in surface temperature as soon as the deconfinement quark mat- 
ter appearing during spin-down. Intermediate magnetic field (10 10 , 10 11 G) 
lead to slight changes in the temperature, but high magnetic field form only 
the temperature plateau at a time. 

In Fig. 6, we present the cooling behavior of different masses stars for mag- 
netic field B=10 12 , 10 n G (left panel) and magnetic field B=10 9 , 10 8 G(right 
panel) with deconfinement heating. The observational data, taken from ta- 
bles 1 and 2 in Page et al.|26], have been shown in left panel. Comparing 



with previous investigation(Kang & Zheng [3j), we find the thermal evolu 
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tion curves of our present work are more compatible with the observational 
data(left panel). In our previous work it seemed that, the NDU processes 
should be triggered easily in the model of relativistic mean field (the criti- 
cal mass for fast cooling occuring is being low )(Glendenning [lij]). In our 
present study, using the APR EOS, NDU processes can not be triggered eas- 
ily in stars which lead to the higher temperatures of the evolution curves 
than in the previous cases. For example, the NDU reactions only appear 
above 1.56 M Q in present model. In the cases of weak magnetic field, stars 
have high temperatures(> 10 5 K) at older ages (> 10 9 yrs). We thus think 
that high temperature of some millisecond pulsars with low magnetic fields 



(Kargaltsev et al [27j]), especially for PSR J0437-4715, can be explained using 
the deconfinement heating model of neutron stars. We can observe that 1.5 
Mq neutron stars follow a similar thermal evolution track as 1.6 M , but 
there is not a period increasing in temperature for 1.7 M Q . Through com- 
paring Fig.6 with Fig. 3, we find the quark matter to appear at the birth of 
star for 1.7M . For 1.5 M & and 1.6 M stars, quark deconfinement occurs 
when the central density gradually increases during spin-down, which results 
in the temperatures of the stars to increase rapidly. This is a characteristic 
signal as quark matter arises during the rotational spin-down of stars for 
weak magnetic case. 

5. Conclusions and discussions 

The chief aim of the present work is to explore the signal of quark matter 
appearing through theoretical simulation of the thermal evolution curves of 
neutron stars with deconfinement heating. We have constructed models of 
rotating neutron stars that follow the mixed phase investigation of Glenden- 
ning based on Hartle's perturbative approach. The total thermal luminosities 
have been obtained using the parameterized approach. 



Recently, Stejner et al [11| have investigated the signature of deconfine- 
ment with spin-down compression in cooling neutron stars. A period of in- 
creasing surface temperature can be produced with the introduction of a pure 
quark core for strongly superfluid stars of strong and intermediate magnetic 
field strength and the latent heat of deconfinement reinforces the signature 
only but is itself relatively less significant. Contrary to their studies, our 
results show that deconfinement heating can drastically affect the thermal 
evolution of neutron stars. The rise of surface temperature of cooling stars, 
as a signature of quark deconfinement, is derived from the deconfinement 
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heating. It is noteworthy that a significant rise of the temperature accompa- 
nies the appearance of quark matter at older ages for low magnetic field stars. 
This may be a evidence for existence of quark matter, if a period of rapid 
heating is observed for a very old pulsar. Deconfinement heating provides a 
new way to study the signal of deconfinement. 

We found that the deconfined signal appears for neutron stars of mass 
1.4M < M < 1.64M Q . The influence of different EOS of the hadron phase 
and the model parameters of the quark phase (bag constant B, coupling 
constant g) on the phase transition densities, rotational structure of neutron 
stars and corresponding internal structure etc have been studied by many 



investigators(Schertler et al.|28j; Pan et al. [29|] ) . The deconfinement heating 
rate and mass range of a deconfinement sinal emerging can be changed with 
varying of these parameters. In future, we will systematically investigate 
the effect of these parameters which which will be the subject of our future 
investigations. 
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Figure 2: Gravitational mass M in solar masses as a function of the central density for 
rotating neutron star configurations with a deconfmement phase transition. The lower 
curve correspond to static configurations, the upper one to those with maximum rotation 
frequency Vk- The lines between both extremal cases connect configurations with the same 
total baryon number. The dotted lines indicate that the quark matter is produced in the 
core of neutron stars. The bag constant of the quark matter is B=108 MeVfm~ 3 
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Figure 3: Central density as a function of rotational frequency for rotating neutron stars 
of different gravitational mass at zero spin. All sequences are with constant total baryon 
number. Dotted horizontal lines indicate that deconfined quark matter is produced and 
dashed horizontal lines indicate that the nucleon direct Urea process is triggered. 
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Figure 4: Thermal evolution curves of a 1.6 M@ neutron star with deconfmement heating 
for various magnetic field strengths (solid curves) and the curves without dcconfincmcnt 
heating (dotted curves) 
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Figure 5: Surface temperature change of a 1.6 M Q neutron star with rotational frequency 
for various magnetic field strengths with deconfinement heating (solid curves) and the 
curves without deconfinement heating (dotted curves) 



15 






1 


















- —-- 1-6M \ 








1-7M S \ 








» PSR J0437-4715 \ 

















Log(t/yr) 



Log(t/yr) 



Figure 6: Thermal evolution curves of neutron stars with deconfinement heating for differ- 
ent stars masses and B=10 12 ,10 11 G(left panel) and B=10 9 ,10 8 G(right panel. Rectangles 
in the left-hand panel indicate observational data on cooling neutron stars with strong 
magnetic fields.) 
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Abstract 

As neutron stars spin-down and contract, the deconfinement phase transi- 
tion can continue to occur, resulting in energy release (so-called deconfinement 
heating) in case of the first-order phase transition. The thermal evolution 
of neutron stars is investigated to combine phase transition and the related 
energy release self-consistently. We find that the appearance of deconfine- 
ment heating during spin-down result in not only the cooling delay but also 
the increase of surface temperature of stars. For stars characterized by in- 
termediate and weak magnetic field strength, a period of increasing surface 
temperature could exist. Especially, a sharp jump in surface temperature can 
be produced as soon as quark matter appears in the core of stars with a weak 
magnetic field. We think that this may serve as evidence for the existence of 
deconfinement quark matter. The results show that deconfinement heating 
facilitates the emergence of such characteristic signature during the thermal 
evolution process of neutron stars. 
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1. Introduction 

Fundamental properties of supranuclear matter in the cores of neutron 
stars, such as the chemical composition and the equation of state, are still 
poorly known. Simulations of the thermal evolution of neutron stars con- 
fronted with soft X-ray, extreme UV, and optical observations of thermal 
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photon flux emitted from their surface provide the most valuable informa- 
tion about the dense matter in the interior of these stars. 

As the interior density gradually increases of neutron stars, deconfine- 
ment phase transition continuously takes place inducing not only structural 
changes but also energy release in case of a first-order phase transition. The 
generation of energy increases the internal energy of the star which is called 
deconfinement heating (Haensel & Zdunik [H,Yu & Zheng (2|,Kang & Zheng 
(if). The thermal evolution of neutron stars is connected with their spin- 
down and the resulting changes in structure and chemical composition (from 
nucleon matter to deconfined quark matter) have been investigated in our 
work. We have investigated the thermal evolution of neutron stars with 
such a deconfinement phase transition(Kang & Zheng (if). The results show 
that deconfinement heating delays dramatically the cooling of neutron stars, 
which have a higher surface temperature compared with traditional cooling 
for the same age stars. 

Many efforts are devoted to explore the observational signal of a decon- 
finement phase transition which have been suggested in the form of char- 
acteristic changes of observables, such as the pulse timing(Glendenning et 
al. j3],Chubarian et al. ji[ ,Poghosyan et al.0]), brightness (Dar &DeRujula 
7|)and surface temperaturefSchaab et al.jsj], Blaschke et al.[9],Yuan & Zhang 
lOj] , Stejner et al. [ll|) of the pulsars during their evolution. The surface 
temperature changes when the quark matter appears in the cores of neutron 
stars. We explore the deconfinement signature by studying the changes of 
surface temperature in the thermal evolution process of neutron stars. 

In this paper, we reinvestigate the thermal evolution of neutron stars and 
look for a characteristic change in the surface temperature with deconfine- 
ment heating. The released energy can also be estimated as a function of 
the change rate of the deconfinement baryon number using the parameterized 
approach. The neutron stars containing quark matter are called hybrid stars. 
We take Glendenning's hybrid stars model (Glendenning (lij ) based on the 
perturbation theory developed by Hartlejl3| to study the rotational evolu- 
tion structure of stars. The Argonne V18 + 5v + UIX* model(APR)(Akmal 



et al. 14|) of hadronic matter and the MIT bag model of quark matter are 



used to construct the model of stars, but the medium effect of quark matter 
has been considered in quasi-particle descript ion (S chert le et al. 15| ) . 



2 



2. Deconfinement phase transition and neutron stars structure 



Early works on deconfinement phase transition, the Maxwell construc- 
tion, show a sharp transition taking place between the two charge-neutral 
hadron and the quark phases(Baym & Chin [la]). In the 1990s, Glenden- 
niiiK (Glendenning 0, H) 

pointed out that this assumption was too re- 



strictive. More generally, the transition can occur through the formation of 
a mixed phase of hadron matter and quark matter, with the total charge 
neutrality being achieved by a positively charged amount of hadronic matter 
and a negatively charged amount of quark matter. Following Glendenning's 
model, we use a standard two-phase description of the equation of state(EOS) 
through which the hadron and quark phases are modelled separately. The 
resulting EOS of the mixed phase is obtained by imposing Gibbs's conditions 
for phase equilibrium with the constraint that the baryon number and the 
electric charge of the system are conserved to the neutron star matter. 

The Gibbs condition for chemical and mechanical equilibrium at zero 
temperature between the two phases reads 

PHp{Pn,Pe) = PQp{Pn,Pe) (1) 

where pnp is the pressure of confined hadron phase and pqp is the pressure 
of deconfined quark phase. 

The conservation laws can be imposed by introducing the quark fraction 
X defined as x — Vq/V. Only two independent chemical potentials remain 
according to the corresponding two conserved charges of the /3-equilibrium 
system. The total baryon number density Pb is 

N B 

PB = ~^ = XPQP + (1 - X)PHP (2) 

the total electrical charge is 



= ^ = XQqp + (1 - x)Qhp (3) 
and the total energy density is 

e = y = X£qp + (1 - x)tHP (4) 

Using the Eqs.(l), (2), (3)and (4), we can obtain the EOS of mixed phase 
matter. In describing the hadronic part of the neutron star, we adopt the 
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APR model (Akmal et al. [14] ). For the EOS, it is based on the models 
for the nucleon interaction with the inclusion of a parameterized three-body 
force and relativistic boost corrections. We use the EOS of an effective mass 



bag-model for the quark matter part of the neutron star(Schertle et al. 15]). 

In Fig.l, we show the model EOS with deconfmement transition, which is 
the typical scheme of a first-order transition. The phase transition construc- 
tion in a two-component system leads to a continuous increase in energy 
per baryon in the mixed phase with increasing density. It is well known 
that hadron matter is the most stable phase at lower densities, and that 
quark matter is the most favorite phase at higher densities. Meanwhile, the 
mixed phase has the lowest energy at intermediate densities. We choose 
the parameters for quark matter EOS with s quark mass m s = 150MeV, 
coupling constant g = 3 and different bag constant B = 85MeVfm~ 3 , 
B = l08MeVfm~ 3 , B = 136MeV 'fm~ 3 respectively. 

With the EOS presented above, we are ready to study the structural 
evolution of the rotating neutron stars. In this paper, we apply Hartle's 
approach (Hartle [l3[) as in Kang & Zheng [3] to investigate the structure 
of the stars. By treating a rotating star as a perturbation on a non-rotating 
star and by expanding the metric of an axially symmetric rotating star in 
even powers of the angular velocity Q, we can obtain the structure of the 
rotating stars. 

The resulting gravitational masses of neutron stars for B = l08MeVfm~ 3 
are shown in Fig. 2, as a function of central baryon density for static stars as 
well as for stars rotating with the maximum rotation frequency v^. The solid 
almost horizontal connect configurations with the same total baryon number. 
In order to explore the increase in central density due to spin down, we cre- 
ated sequences of neutron star models. A Model in particular sequence has 
the same constant baryon number, increasing central density and decreasing 
angular velocity. 

3. Deconfmement heating 

There is deconfmement heating production due to spin-down in neutron 
stars due to the nuclear matter continuously converting into quark matter. 
The released energy had been estimated as a function the rate of change 
the deconfmement baryon number using the parameterized approach(Kang 
& Zheng [3|). Recently we studied the mechanism of energy release in de- 



tail (Kang et al. [18] ). Through studying a random process of infinitesimal 
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compression for the mixed phase region, we can calculate the energy release 
per baryon using the following formula 



Se - 5e = (-p- - -^—)5p B , (5) 

OpB 0p B 

where denotes the enthalpy change per baryon . 

The deconfmement heating is coupled with the rotation evolution of neu- 
tron stars. Combining the energy change with the evolutionary structure of 



neutron stars, we get the total heat luminosity (Kang et al. 18) 19) ) 



f de 

H dec = I -^v(t)p B dV (6) 

where v is the rotation frequency of the star. The spin-down of stars is due 
to the magnetic dipole radiation. The rotation frequency is given by 

v = -^L^Vsin 2 ^ (7) 
31c 6 

where / is the stellar moment of inertia, p = \BB? is the magnetic dipole 
moment, and 6 is the inclination angle between the magnetic and rotational 
axes. We now combine the energy release behavior with spin-down. Our 
recent work shows that exact calculations agree well with the earlier order 
of magnitude estimates and supports the previous parameterized approach. 
Because the change of rotating stars is sufficiently slow, the formula of heat 
luminosity can be replaced by the simple parameterized form(Kang & Zheng 

!), 

H dec = q^vit) (8) 
dv 

where q, about O.lMeV, is the mean value of energy release in the mixed 
phase and Nq represents the baryon number of quarks in the interior of the 
star. 



4. Signal of quark deconfinement and thermal evolution of neutron 
stars 

The cooling of neutron stars could take place via two channels - neutrino 
emission from the entire star and thermal emission of photons through the 
transport of heat from the internal layers to the surface. Neutrino emission 
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is generated in numerous reactions in the interior of neutron stars, e.g. as 
reviewed by Page et al.j2o[. For the calculation of cooling of the hadronic 
part of the neutron star, we use the main processes including the nucleon 
direct Urea (NDU) and the nucleon modified Urca(NMU) and the nucleon 
bremsstrahlung(NB). For the quark matter, we consider the main process: 
the quark direct Urea (QDU) processes on unpaired quarks, the quark modi- 
fied Urea (QMU) and the quark bremsstrahlung(QB). For pure neutron stars, 
the direct Urea reaction (the most efficient) is allowed only at very high den- 
sities because it is impossible to satisfy the conservation of momentum unless 
the proton fraction exceed the value where both the charge neutrality and 
the triangle inequality can be observed (Lattimer et al. [21] ) . However, for 
neutron stars that contains quark matter, this is not so because charge neu- 
trality does not have to be conserved locally for a mixed phase. Hence the 
NDU process is active in the mixed phase. 

The calculation of the evolution of the thermal energy of neutron stars 
(heating and cooling) is achieved by coupling their rotating structure and de- 
confinement phase transition. Fig. 3 displays the central density of different 
masses neutron stars for B = 108MeVfm~ 3 , as a function of their rotational 
frequency. In Fig.3, the dotted horizontal line indicates the deconfined quark 
matter produced and dashed horizontal line indicate the NDU processes trig- 
gered. In the interior of these stars, the emerging of deconfined quark matter 
accompanies the gradual energy release leads to rising of surface tempera- 
ture, and appearance of the NDU process can result in a rapid decrease of 
the temperature during the spin-down. For M = 1.5, 1.55, 1.6M Q neutron 
stars, the deconfinement phase transition occurs during the thermal evolu- 
tion process which may lead to appearance of a characteristic signature. We 
will now discuss in detail the quark deconfinement signal. 

We combine the equation of thermal balance with the rotating structure 
equations of the stars(Kang & Zheng j3|, Hartle [l3j]) and rewrite the energy 
equation in the approximation of an isothermal interior(Glen Sz Sutherland 

0) 

CviTuv)^ = -L?{Ti,v) - L~(T.,v) (9) 

C v (T u v) = c(r,T)(l - ^M)-V2 47rr 2 dr (1Q) 

Jo r 

L?(T h v)= j R{V \{r,T){l- 2 -^)-y^H^dr (11) 
Jo r 
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Where T s is the effective surface temperature, Tj(£) = T(r, t) is the redshifted 
internal temperature; T(r, t) is the local internal temperature of matter, 
and $(r) is the metric function(describing gravitational redshift)(Yakovlev & 
Haensel |23[). Furthermore, L£°(Tj, t>)and CV(Tj,t>) are the total redshifted 
neutrino luminosity and the total stellar heat capacity, respectively, which 
are functions of rotation frequency and temperature; c(r, T) is the heat ca- 
pacity per unit volume. = AnR 2 {v)aTf{l — R g /R) is the surface photon 
luminosity as detected by a distant observer(i? ff is the stellar gravitational 
radius). The effective surface temperature which is detected by a distant 
observer is T s °° = T s a/I — R g /R. T s is obtained from the internal temper- 



ature by assuming an envelope model(Gudmundsson et al. [24|, Potekhin 
et al |25|). Using to Eqs.(7)-(11), we can simulate the thermal evolution of 
neutron stars with deconfinement heating. 

In Fig.4 and Fig.5, we present thermal evolution behavior of a 1.6M Q 
neutron star for different magnetic fields (10 9 — 10 12 G). Due to the coupling 
of thermal evolution and spin-down, all curves(with deconfinement heating 
and without deconfinement heating) show clear magnetic field dependence. 
During the star's spin down, deconfined quark matter appears in the core of 
the star at a spin frequency of v — 1123Hz, the surface temperature drops 
rapidly when the NDU process(enhanced cooling) occurs at a spin frequency 
of v = 492Hz. It is evident that the temperature of the curves with de- 
confinement heating(solid curves) are higher than for the standard cooling 
scenario without deconfinement (dotted curves). We can observe a competi- 
tion between cooling and heating processes from the heating curves, where 
deconfinement heating can produce a characteristic rise of surface tempera- 
ture and even dominate the history of thermal evolution. Eventually, they 
reach a thermal equilibrium, where the heat generated is radiated away at 
the same rate from the star surface. We find the weaker magnetic fields have 
the larger change of temperature. The low magnetic field (10 9 G) produces a 
sharp jump in surface temperature as soon as the deconfinement quark mat- 
ter appearing during spin-down. Intermediate magnetic field (10 10 , 10 11 G) 
lead to slight changes in the temperature, but high magnetic field form only 
the temperature plateau at a time. 

In Fig. 6, we present the cooling behavior of different masses stars for mag- 
netic field B=10 12 , 10 n G (left panel) and magnetic field B=10 9 , 10 8 G(right 
panel) with deconfinement heating. The observational data, taken from ta- 
bles 1 and 2 in Page et al.|26], have been shown in left panel. Comparing 



with previous investigation(Kang & Zheng [3j), we find the thermal evolu 



7 



tion curves of our present work are more compatible with the observational 
data(left panel). In our previous work it seemed that, the NDU processes 
should be triggered easily in the model of relativistic mean field (the criti- 
cal mass for fast cooling occuring is being low )(Glendenning [lij]). In our 
present study, using the APR EOS, NDU processes can not be triggered eas- 
ily in stars which lead to the higher temperatures of the evolution curves 
than in the previous cases. For example, the NDU reactions only appear 
above 1.56 M Q in present model. In the cases of weak magnetic field, stars 
have high temperatures(> 10 5 K) at older ages (> 10 9 yrs). We thus think 
that high temperature of some millisecond pulsars with low magnetic fields 



(Kargaltsev et al [27j]), especially for PSR J0437-4715, can be explained using 
the deconfinement heating model of neutron stars. We can observe that 1.5 
Mq neutron stars follow a similar thermal evolution track as 1.6 M , but 
there is not a period increasing in temperature for 1.7 M Q . Through com- 
paring Fig.6 with Fig. 3, we find the quark matter to appear at the birth of 
star for 1.7M . For 1.5 M & and 1.6 M stars, quark deconfinement occurs 
when the central density gradually increases during spin-down, which results 
in the temperatures of the stars to increase rapidly. This is a characteristic 
signal as quark matter arises during the rotational spin-down of stars for 
weak magnetic case. 

5. Conclusions and discussions 

The chief aim of the present work is to explore the signal of quark matter 
appearing through theoretical simulation of the thermal evolution curves of 
neutron stars with deconfinement heating. We have constructed models of 
rotating neutron stars that follow the mixed phase investigation of Glenden- 
ning based on Hartle's perturbative approach. The total thermal luminosities 
have been obtained using the parameterized approach. 



Recently, Stejner et al [11| have investigated the signature of deconfine- 
ment with spin-down compression in cooling neutron stars. A period of in- 
creasing surface temperature can be produced with the introduction of a pure 
quark core for strongly superfluid stars of strong and intermediate magnetic 
field strength and the latent heat of deconfinement reinforces the signature 
only but is itself relatively less significant. Contrary to their studies, our 
results show that deconfinement heating can drastically affect the thermal 
evolution of neutron stars. The rise of surface temperature of cooling stars, 
as a signature of quark deconfinement, is derived from the deconfinement 
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heating. It is noteworthy that a significant rise of the temperature accompa- 
nies the appearance of quark matter at older ages for low magnetic field stars. 
This may be a evidence for existence of quark matter, if a period of rapid 
heating is observed for a very old pulsar. Deconfinement heating provides a 
new way to study the signal of deconfinement. 

We found that the deconfined signal appears for neutron stars of mass 
1.4M < M < 1.64M Q . The influence of different EOS of the hadron phase 
and the model parameters of the quark phase (bag constant B, coupling 
constant g) on the phase transition densities, rotational structure of neutron 
stars and corresponding internal structure etc have been studied by many 



investigators(Schertler et al.|28j; Pan et al. [29|] ) . The deconfinement heating 
rate and mass range of a deconfinement sinal emerging can be changed with 
varying of these parameters. In future, we will systematically investigate 
the effect of these parameters which which will be the subject of our future 
investigations. 
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Figure 2: Gravitational mass M in solar masses as a function of the central density for 
rotating neutron star configurations with a deconfmement phase transition. The lower 
curve correspond to static configurations, the upper one to those with maximum rotation 
frequency Vk- The lines between both extremal cases connect configurations with the same 
total baryon number. The dotted lines indicate that the quark matter is produced in the 
core of neutron stars. The bag constant of the quark matter is B=108 MeVfm~ 3 
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Figure 3: Central density as a function of rotational frequency for rotating neutron stars 
of different gravitational mass at zero spin. All sequences are with constant total baryon 
number. Dotted horizontal lines indicate that deconfined quark matter is produced and 
dashed horizontal lines indicate that the nucleon direct Urea process is triggered. 
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Figure 4: Thermal evolution curves of a 1.6 M@ neutron star with deconfmement heating 
for various magnetic field strengths (solid curves) and the curves without dcconfincmcnt 
heating (dotted curves) 
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Figure 5: Surface temperature change of a 1.6 M Q neutron star with rotational frequency 
for various magnetic field strengths with deconfinement heating (solid curves) and the 
curves without deconfinement heating (dotted curves) 
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Figure 6: Thermal evolution curves of neutron stars with deconfinement heating for differ- 
ent stars masses and B=10 12 ,10 11 G(left panel) and B=10 9 ,10 8 G(right panel. Rectangles 
in the left-hand panel indicate observational data on cooling neutron stars with strong 
magnetic fields.) 
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